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Peristaltic flow of a Hyperbolic tangent fluid through 
a porous medium in a vertical channel 
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Abstract — In this paper, we investigated the peristaltic 
transport of a hyperbolic tangent fluid through a porous 
medium in a vertical channel under the assumption of long 
wavelength. The expression for the velocity and axial pressure 
gradient are obtained by employing perturbation technique. The 
effects of various pertinent parameters on the time-averaged 
flow rate are discussed with the help of graphs. 

Index Terms — Darcy number, Hyperbolic tangent fluid, 
peristalsis. Vertical Channel. 


I. Introduction 

A Peristaltic pump is a device for pumping fluids, generally 
from a region of lower to higher pressure, by means of a 
contraction wave traveling along a tube-like structure. This 
traveling wave phenomenon is referred to as ‘Peristalsis’. 
Peristalsis originated naturally as a means of pumping 
physiological fluids from one place in the body to another, 
and is the primary pumping physiological fluids from one 
place in the body to another, and is the primary pumping 
mechanism in swallowing (and indeed all the way through the 
alimentary canal) in the ureter, the bile duets, the ductus 
efferentes of the male reproductive tract, and even in some 
small blood vessels. Humankind has borrowed the idea and 
used it in applications where the material being pumped must 
not be contaminated (eg blood) or is corrosive and should not 
be in contact with the moving parts of ordinary pumping 
machinery. Analytical solutions were obtained for peristaltic 
flows by assuming either small amplitude but arbitrary 
Reynolds number [1] -[2] or arbitrary amplitude with small 
curvature and negligible inertia [3]. 

It is well known that some fluids which are encountered in 
chemical applications do not adhere to the classical 
Newtonian viscosity prescription and are accordingly known 
as non-Newtonian fluids. One especial class of fluids which 
are of considerable practical importance is that in which the 
viscosity depends on the shear stress or on the flow rate. The 
viscosity of most non-Newtonian fluids, such as polymers, is 
usually a nonlinear decreasing function of the generalized 
shear rate. This is known as shear-thinning behavior. Such 
fluid is a hyperbolic tangent fluid [4]. Reference [5] have first 
investigated the peristaltic flow of a hyperbolic tangent fluid 
in an asymmetric channel. Reference [6] have analyzed the 
peristaltic transport of a Tangent hyperbolic fluid in an 
endoscope numerically. Reference [7] have discussed the 
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peristaltic flow of a hyperbolic tangent fluid in an inclined 
asymmetric channel with slip and heat transfer. 

Flows through a porous medium occur in filtration of fluids. 
Hall effects on peristaltic flow of a Maxwell fluid in a porous 
medium were investigated by [8]. Reference [9] have 
discussed the effect of magnetic field on the peristaltic motion 
of a Carreau fluid through a porous medium with heat 
transfer. Peristaltic motion of a couple stress fluid through a 
porous medium in a channel with slip condition was studied 
by [10]. Reference [11] have studied the peristaltic MHD 
flow of a Bingham fluid through a porous medium in a 
channel. Peristaltic flow of a non-Newtonian fluid through a 
porous medium in a tube with variable viscosity using 
Adomian decomposition method was investigated by [12]. 

In view of these, we studied the peristaltic flow of a 
hyperbolic tangent fluid through a porous medium in an 
inclined channel under the assumption of long wavelength. 
The expression for the velocity and axial pressure gradient are 
obtained by employing perturbation technique. The effects of 
various pertinent parameters on the time-averaged flow rate 
are discussed with the help of graphs. 

II. Formulation of the Problem 

We consider the peristaltic motion of a hyperbolic tangent 
fluid through a porous medium in a two-dimensional 
symmetric vertical channel of width 2 CL . The flow is 
generated by sinusoidal wave trains propagating with constant 
speed C along the channel walls. Fig. 1 illustrates the 
schematic diagram of the channel. 

The wall deformation is given by 

2 77- 

Y = ±H(X,t) = ±a±b cos — ( X-ct ), (2.1) 

X 

where b is the amplitude of the wave, A - the wave length and 
X and Y - the rectangular co-ordinates with X measured 
along the axis of the channel and Y perpendicular to X . Let 
(U,V) be the velocity components in fixed frame of 
reference (X, Y ) . 

The flow is unsteady in the laboratory frame (A,T) . 
However, in a co-ordinate system moving with the 
propagation velocity c (wave frame (x, y)), the boundary 
shape is stationary. The transformation from fixed frame to 
wave frame is given by 

x = X —ct,y — Y,u — U — c,v — V (2.2) 

where ( u,v ) and ( U,V ) are velocity components in the wave 
and laboratory frames respectively. 
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Fig. 1 The physical model 
The constitutive equation for a Hyperbolic Tangent 

fluid is 

r = -[*7« 0 +(%+/ 7 „)tanli( r j>)'’]x (2.3) 

where T is the extra stress tensor, Tj^ is the infinite shear 
rate viscosity, TJ 0 is the zero shear rate viscosity, F is the 
time constant, fl is the power-law index and y is defined as 
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where p is the density and k is the permeability of the 
porous medium. 

The corresponding dimensional boundary conditions are 
u = -c at y = H (2.9) 


du 

— = 0 at y = 0 

dy 


( 2 . 10 ) 

Introducing the non-dimensional variables defined 


by 

_ x _ y _ u _ v „ a 

x = — , y = — , w = — , v = — , o = — , 

X a c cS A 

_ pa 2 . b , H - ct 

P = — ,</> = — ,h = — , t = — , 

77 0 c/l a a A 
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into the Equations (2.6) - (2.8), reduce to (after dropping the 
bars) 
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where n is the second invariant stress tensor. We consider in 
the constitutive equation (2.3) the case for which ij^- 0 and 
r> <1, so the Eq. (2.3) can be written as 

z ‘ = -%( r b)V = -%( 1+r r-i)V = -%( 1 +«[ r r-i])r 

(2.5) 

The above model reduces to Newtonian for T = 0 and 
n = 0 . 

The equations governing the flow in the wave frame of 
reference are 


Re is the 


k . 


Reynolds number, Fr is the Froude number and Da - — is 

a 

the Darcy number. 

Under long wavelength approximation (£ — » 0) , the Eqs. 
(2.13) and (2.14) become 

1 Rp 

— (k + 1) + — (2.15) 
Da V J Fr 

(2.16) 


dp d J 

l + n 

( Mr du y 

du\ 

dx dy | 


l dy J. 

Fl 


^ = 0 

dy 


From Eq. (2.15) and (2.16), we get 

, 2 . 


dp , sdu d 

— = ( 1 - n ) — - + n We — 
dx v J dy 2 dy 


r du V 
\dy j 


1 / ^ Re 
— (u- bl)-b — 
Da v ’ Fr 


(2.17) 

The corresponding non-dimensional slip boundary 
conditions in the wave frame are given by 
u — —\ at y = h = l + ^cos 2 ^x (2.18) 

du 


dy 


= 0 


at y = 0 


(2.19) 
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The volume flow rate q in a wave frame of 
reference is given by 

h 

q = J udy . (2.20) 

o 

The instantaneous flow Q(X,t) in the laboratory frame is 

h h 

Q(X,t) = j UdY = | (w + i)dy = q+h (2.21) 

0 0 

The time averaged volume flow rate Q over one period 

f >0 

T = — of the peristaltic wave is given by 

III. Solution 

Since Eq. (2.17) is a non-linear differential equation, so it is 
not possible to obtain closed form solution. Therefore we 
employ regular perturbation to find the solution. 

dp 

For perturbation solution, we expand w,— and q as 

dx 


1 ( 

dPo 

Re'] 

sinh ah - ah cosh ah 

<r 3 (l-n)l 

v dx 

Fr J 

cosh ah 


% 

(3.11) 

From Eq. (3.1 1), we have 

dp 0 (q 0 + h)a 3 (l-ft) cosh ah 
dx 


Re 


(3.12) 


[ sinh crh - crh cosh crh] Fr 

3.4 Solution for system of order We 1 

Substituting Eq. (3.10) in the Eq. (3.7) and solving the Eq. 
(3.7), using the boundary conditions (3.8) and (3.9), we obtain 


J* Qdt — q + 1 

(2.22) 

u, = 


f dp^_Re Y 

0 



n 

^ dx Fr) 


follows 

u = u 0 +Weu,+ O (We 2 ) 


— = + We — + O (We 2 ) 

dx dx dx v 7 

q = q 0 +Weq l +O(We 2 ) 


(3.1) 

(3.2) 

(3.3) 


Substituting these equations into the Eqs. (2.17) - (2.19), 
we obtain 

3.1. System of order We° 


dp ( 
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and the respective boundary conditions are 
u () = - 1 at y-h 

du. 


(3.5) 
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ay 
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and the respective boundary conditions are 
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3.3 Solution for system of order We° 

Solving Eq. (3.4) using the boundary conditions (3.5) and 
(3.6), we obtain 
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Uc =- 
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Re 
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-1 


a 2 (l -n) , 
where cr = 1/ ^Da(l-n) . 

The volume flow rate q 0 is given by 


(3.10) 
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The volume flow rate q, is given by 
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where 

A = 4- cosh crh + 2 sinh Icih sinh crh — cosh ah cosh 2a h . 
From Eq. (3.14) and (3.12), we have 
q,a 3 (l- ft) cosh ah 


dPi 
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and using the relation = — ~We -FL an( i neglecting 


Substituting Equations (3.12) and (3.15) into the Eq. (3.2) 

^-We^ 

dx dx dx 
terms greater than 0(We ) , we get 

^ (q + h)a 3 (l-ft) cosh ah ^ 


dp 

dx 


[sinh ah - ah cosh ah] 
We Ana 5 (# + ^) 


Re 


(3.16) 


v 6 (sinh ah -ah cosh ah) Fr 

The dimensionless pressure rise per one wavelength in the 
wave frame is defined as 

A p=\'^-dx (3.17) 

io dx 

Note that, as Da —> oo , We — » 0, Re — > 0 and n — > 0 our 
results coincide with the results of Shapiro et al. (1969). 

IV. Discussion of the Results 
Fig. 2 illustrates the variation of pressure rise A p 

with time-averaged volume flow rate Q for different values 
of We with <f> = 0.5 , ft = 0.5 , Re = 1 , Fr = 0.2 and 
Da = 0.01 . It is found that, the time- averaged volume flow 
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rate Q increases with increasing We in all the three regions 
(namely, pumping (Ap > 0) , free-pumping (Ap = 0) and 
co-pumping (Ap<0) regions). 

The variation of pressure rise Ap with 

time-averaged volume flow rate Q for different values 
of n with (j> = 0.5 , We = 0.05 , Re = 1 , Fr = 0.2 and 
Da = 0.01 is shown in Fig. 3. It is observed that, the 
time-averaged volume flow rate Q decreases with an 
increase in n in both the pumping and free pumping regions, 
while it increases with increasing n in co-pumping region for 
chosen Ap(<0). 

Fig. 4 depicts the variation of pressure rise Ap with 

time-averaged volume flow rate Q for different values of 
Da with (j) = 0.5 , n = 0.5 , Re = 1 , Fr = 0.2 and 
We = 0.05 . It is noted that, the time-averaged volume flow 
rate Q decreases with increasing Da in the pumping region, 
while it increases with increasing Da in both the 
free-pumping and co-pumping regions. 

The variation of pressure rise Ap with 

time-averaged volume flow rate Q for different values 
of ([> with We = 0.05 , n = 0.5 , Re = l , Fr = 02 and 
£><2 = 0.01 is depicted in Fig. 5. It is found that, the 
time-averaged volume flow rate Q increases with an 
increase in (j) in both the pumping and free-pumping regions, 
while it decreases with increasing $ in the co-pumping 
region for chosen Ap(< 0) . 

Fig. 6 illstrates the variation of pressure rise Ap 

with time-averaged volume flow rate Q for different values 
of Re with <f> = 0.5 , n = 0.5 , We = 0.05 , Fr = 0.2 and 
Da = 0.01 . It is observed that, the time-averaged volume 
flow rate Q increases with increasing Reynolds number Re 
in all the three regions. 

The variation of pressure rise Ap with 

time-averaged volume flow rate Q for different values of 
Fr with (f) = 0.5 , n = 0.5 , Re = 1 , W<? = 0.05 and 
£><2 = 0.01 is shown in Fig. 7. It is found that, the 
time-averaged volume flow rate Q decreases with increasing 
Froud number Fr in all the three regions. 

V. Conclusions 

In this paper, we studied the peristaltic flow of a hyperbolic 
tangent fluid through a porous medium in a vertical channel 
under the assumption of long wavelength. The expression for 
the velocity and axial pressure gradient are obtained by using 
perturbation technique. It is found that in the pumping region 

the time-averaged flow rate Q increases with increasing We , 
(j) and Re , while it decreases with increasing n , Da and 
Fr . Further, it is observed that the pumping is more for 
hyperbolic tangent fluid than that of Newtonian fluid. 



Q 


Fig. 2 The variation of pressure rise Ap with time-averaged 

volume flow rate Q for different values of We with 
^ = 0.5, ?7 = 0.5, Re = 1 , Fr = 0.2 and £><2 = 0.01. 



Fig. 3 The variation of pressure rise Ap with time-averaged 

volume flow rate Q for different values of n with (f) = 0.5 , 
We = 0.05 , Re = 1 , Fr = 0.2 and Da = 0.01 . 



Fig. 4 The variation of pressure rise Ap with time-averaged 

volume flow rate Q for different values of Da with 
(/) = 0.5 , n = 0.5 , Re = 1 , Fr = 0.2 and We = 0.05 . 
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Fig. 5 The variation of pressure rise A p with time-averaged 

volume flow rate Q for different values of (j) with 
We = 0.05 , n = 0.5 , Re = 1 , Fr- 0.2 and Da = 0.01 . 



Q 


Fig. 6 The variation of pressure rise A p with time-averaged 

volume flow rate Q for different values of Re with 
^ = 0.5, rc = 0.5, We = 0.05 , Fr = 0.2 and Da = 0.01. 



Q 

Fig. 7 The variation of pressure rise A p with time-averaged 

volume flow rate Q for different values of Fr with (/) = 0.5 , 
n = 0.5 , Re = 1 , We = 0.05 and Da = 0.01 . 
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